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ABSTRACT 

We investigate the relation between the regularized multi-propagator method, called "Reg 
PT", and the standard perturbation theory. Reg PT is one of the most successful models to 
describe nonlinear evolution of dark matter fluctuations. However, Reg PT is a mathemati- 
cally unproven interpolation formula between the large-scale solution calculated by the standard 
perturbation theory and the limiting solution in the small scale calculated by the multi-point 
propagator method. In this paper, we give an alternative explanation for Reg PT in the context 
of the standard perturbation theory, showing that Reg PT does not ever have more effective in- 
formation on nonlinear matter evolution than the standard perturbation theory. In other words, 
the solutions of the standard perturbation theory reproduce the results of A^-body simulations 
better than those of Reg PT, especially at the high-A: region. This fact means that the standard 
perturbation theory at the two-loop level is still one of the best predictions of the nonlinear power 
spectrum to date. Nevertheless, the standard perturbation theory has not been preferred because 
of the divergent behavior of the solution at small scales. To solve this problem, we also propose 
a modified standard perturbation theory which avoids the divergence. 

Subject headings: dark matter — large-scale structure of Universe 



Introduction 



What is the best model to predict the nonlinear matter power spectrum ? In this paper, we show 
that the standard perturbation (SPT) theory with two-loop corrections is still one of the best models, even 
though various modified perturbation theories have been proposed in the last 5-10 yr. In doing so, we give 
the relation between the SPT and the regularized multi-propagator method. 

One of the key quantities in modern cosmology is the matter power spectrum because it contains a lot 
of important information on evolution of the universe as well as on structure formation. Recent observation 
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of the baryon acoustic oscillations (BAO) in the power spectrum provides a new method to precisely restrict 
cosmological parameters (Eisenstein et al. 2005). 

Recent progress on cosmological observation greatly motivates various theoretical studies of accurate 
calculation of the nonlinear matter power spectrum. At present, one of the most successful models is given 
by the multi-point propagator method (F-expansion method) with the regularized treatment of propagators 
presented by Bernardeau et al. (2008, 2012), called "Reg PT" in Taruya et al. (2012). In particular, Taruya 
et al. (2012) computed the solution of Reg PT with two-loop corrections, giving theoretical predictions of 
the power spectrum which agree well with A^-body simulations at BAO scales. 

Various modified perturbation theories partially sum up high-order terms in the SPT, called "resum- 
mation theories" . Since SPT has exact but formal solutions at any order in perturbation theory, any 
resummation theory should be represented in the context of SPT. The coefficients of the F-expansion can 
easily be represented by the kernel function in SPT. This is why we consider the relation between SPT and 
the F-expansion method. 

We also proposed a method called the "Wiener Hermite (WH) expansion" which gives an accurate power 
spectrum over relatively wide range of wavenumber (Sugiyama & Futamase 2012). There, we established 
the relation between the WH expansion and SPT, and showed the equivalence between the WH expansion 
and the F-expansion method (Bernardeau et al. 2008, 2012; Taruya et al. 2012). Furthermore, by using an 
approximation of the kernel functions in SPT, we derived the exponentially damping behavior of the power 
spectrum at each order of the WH expansion (F-expansion) which has been known in the renormalized 
perturbation theory (RPT; Crocce & Scoccimarro (2006a, b, 2008)). 

We extend the approximation method of the kernel functions introduced in our previous work (Sugiyama 
& Futamase 2012). As a result, we give an alternative explanation for Reg PT, and propose a natural 
extension of Reg PT. The extended Reg PT gives the SPT solution with negligible correction terms. Thus, 
we find that the Reg PT solution does not ever have more effective information on the evolution of dark 
matter than the SPT solution. In addition, we also propose a modified version of the SPT solution which 
avoids the divergent behavior of the SPT solution at small scales. 

The A-body simulation results used in this paper were presented by Valageas & Nishimichi (2011). These 
results and initial conditions at Zini = 99 were created by the public A^-body codes GADGET2 and 2LPT 
code, respectively (Springel 2005; Crocce et al. 2006). These A^-body simulations contain 2048'^ particles 
and were computed by combining the results with different box sizes 2048/i~^ Mpc and 4096 /i~^Mpc, called 
Lll-A^ll and il2-A^ll. The cosmological parameters we used were presented by the Wilkinson Microwave 
Anisotropy Probe five year release (Komatsu et al. (2009), D.^ = 0.279, f^A = 0.721, = 0.046, h = 0.701, 
Us ~ 0.96 and erg = 0.817). We used the program which is available on Taruya's homepage to compute the 
predicted power spectra ^. 

This paper is organized as follows. In Section 2, the relation between SPT and the F-expansion is 
reviewed. Section 3 proves the approximation of the kernel functions. Using the approximated kernel 
functions, the alternative explanation for Reg PT is given in Section 4. In Section 5, it is proposed that 
the natural extension of Reg PT is equivalent to the solution of SPT with some correction terms. We also 
present the modified version of the SPT solution in Section 6. There, we reconfirm that the two-loop SPT 
solutions are better than those of the two-loop Reg PT and the closure theory, and show that the divergent 
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behavior of the solution in SPT is indeed removed by our modification. In Section 7, we summarize and 
discuss our results. 



2. Review of SPT and the F-expansion method 

1 /2 

We shall briefly review SPT. In SPT, we choose a cosmological model with / = rtm where / = 
din D/ din a is the linear growth rate with D and a being the linear growth factor and the scale factor. As 
shown by Lahav et al. (1991), this model is a good approximation in practice. In this situation, the density 
perturbation and the velocity divergence of dark matter are expanded as (Bernardeau et al. 2002), 

oo oo 

d{z, k) = ^ D^S,,{k), 0{z, k) = -aHf ^ i?"0„(k), (2-1) 

n— 1 n— 1 

where H is the Hubble parameter, and (5„ and 9n are time-independent quantities with a Fourier mode k. 
The velocity divergence is defined as = V ■ v, where v is the velocity of dark matter. The nth-order 
solutions in SPT are given by, 

^n(k) = J ^^■■■J ^^(27r)^(5D(k-p[i,„])i^„([pi,p„])(5L(pi)---(5L(p«), 

^"^^^ = |^---y'^(2^)''5D(k-p[i,„])G„([pi,p„])<5L(pi)---<5L(p„), (2-2) 

where P[i.ri] = pi + ■ • • + p,i, Fn{pi, . . . , p„) = -F')i([Pi, Pn]), and (5d is the Dirac delta function. Linearized 
quantities, such as (5l and 6'l, are denoted by the subscript L. The kernel functions F and G are constructed 
from the mode coupling functions a(ki,k2) = (ki -I- k2) • k2/fci and /3(ki,k2) = |ki -I- k2p(ki • k2)/2fc^fcf 
according to the following recursion relation: 

n-l 

nK([pi,p„]) - G„([pi,p„]) = ^ Gm([qi,qm])Q;(q[i 

,m] 1 *l.[?n+l,n] )-F'n-m([qm+l, Qn]), 



m— 1 

3 2n + 1 "^^ 

--^Fni[Pl,Pn]) -l ^ G„([pi,p„]) = ^ Grni[qi,Clm])f3{ci[l ,m\ 7 

m— 1 

(2-3) 

where Fi ^ Gi = 1, and the wavenumber vectors {qi, . . . , qn} = {pi, . • . , Pn} iu the right hand side are 
completely-symmetrized. 

The relation between the rth-order coefficient in the F-expansion and the kernel function is given by 
(Bernardeau et al. 2008; Sugiyama & Futamase 2012) 



Ft--) (Z, [ki, k.]) EE D'r[:li[k,,kr]) + J2 Z?'-+'"ri:Lp([ki, k.]), (2-4) 

n=l 



with 



r!:Lp([ki,k.]) 

_ 1 {r + 2n)l f d^pi f d^pn 



r\ 2"n! 



/d pi f d Pji 

(2^ "J (^:;;^^'-t-2"(^l''^'-]'Pl' ~P1' ■ • ■ ,Pn,-Pn)-PL(Pl) • --PLiPn 

(2-5) 
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whcrc Pl is the present linear power spectrum. The n-loop correction to the rth-order coefficient of the 
F-expansion r|^'^|^^p is of ©((Pl)"), and involves n-dimensional integrals. The full nonlinear matter power 
spectrum is given by 

oo 

P{z,k)^Y.Pr\^:k), (2-6) 

r=l 

where Pp is the rth-order contribution to the power spectrum in the F-expansion, defined as 

p('^)(z,fc)EEH J y ^(2^)35^(k-k[i,.])[FW(z,[ki,k.])] PL(fci)---PL(fc.). (2-7) 

Note that Pr which comprises Fl^^^J^^p x f|j^_\qqp is of O ((Pl)*^^"^™), and has {r + n + m — l)-dimensional 
integrals, also called the {r + + m — l)-loop correction term. 

The velocity divergence satisfies the same relation as the above by replacing F with G. 



3. Approximation of Kernel Functions 

Since the F-expansion method is rewritten by the kernel functions in SPT, the problem to solve the 
nonlinear evolution of dark matter eventually reduces to compute the kernel functions. Therefore, we develop 
an approximation of the kernel functions. 

We shall impose the condition that the amplitude of a Fourier mode p in P„ and Gn is much smaller 
than that of the others, 

|p,|>|p|^0, for {^ = (3-1) 
In this condition, we prove the following approximation of the kernel functions: 

^n([Pl,P«-l],p) ^ - ( ^ ) P„_i([pi,p„_i]), 

G„([Pi,P„-i],p) ^ -f ^'''"','''^ )G„,i([pi,p„_i]). (3-2) 

We provide proof of this by the recursion relation in Equation (2-3) and induction in n. Note that 
these expressions are slightly different from and more generalized than those presented in our previous 
work (Sugiyama & Futamase 2012). 

First, the kernel functions F2 and G2 clearly satisfy Equation (3-2) for |p2| 3> |pi| ^ 0: 
p. n ^1 - 5 I pi - Pa /Pi P2\ 2 (pi -pz)' 1 / Pi ■P2 \ 

G2(P1,P2)| , ^ 3 Ipi-p. M 4 (p,^p,)^ 1 /p^p,\ 

The next step is to prove the following relation for |p| — > 0: 



^ri+l([Pl,Pn],P) ^ _J_ I P[l-"1 P j P„([pi,p„]), 

G„+i([pi,p„],p) ^ ^ ("EM-P )G„([pi,p„]), (3-4) 
n + I \ p-^ ' 
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with the assumption that Fi and Gi for {i ~ 1, . . . ,n} satisfy Equation (3-2). For this purpose, we rewrite 
the first fine in Equation (2-3) as follows: 

(n + l)F„+i([pi,p„],p) - G„+i([pi,p„],p) 



ji-i 



( ^2+1 ) <^™+l([qi'q™]'P)"(q[l,™] +P.q[m+l,n])-P'«-m([qm-M,q«]) 



7n—l 
n-1 



+ ( — — ) ^™([qi'q™])"(q[i,™]'q["i-Hi,"] +p)K-Hi-m([qm-Hi,qn],p) 

H 7T"(P' P[l.n])^»([Pl' P"]) -I l-r"(P[l.n]' P)<^"([Pl' P'l])- 

n+1 nH-1 
This equation is approximated as 

(n-t- l)F„+i([pi,p„],p) - G„+i([pi,p„],p) 



1 / P[l,n] • P 



Gm([qi, q,„])a(q[i^,„],q[,m+i,„])-F'n-m([qm+i, qn]) 



n + 1 \ 

1 ^Pd.npP^ i^„([pi,P„]) + :^G„([Pl,P„] 



n + 1 \ p- 
where we used the following relations for |p| — > 0: 



n+1 



Gm+i([qi,qm],p) 

-Fn+l-m([qm+l, qn], P) 
a(q[i,m] + P, q[m+l,ji]) 



2 — ) Gm([qi,q,i 



m + 1 
1 



q[m+l,n] • P 



a(q[l,m],q[m+l.n] + P) 
"(P[l,n],P) 

Using the first line in Equation (2-3) and ignoring the last term in Equation (3-6), we find 



n + 1 — m \ 

"(q[l,Tn],q[m-|-l,n]), 
a(q[l,m],q[m+l,n]), 
1. 



1 ^n—m ( [Qtti ; Qn] ) : 



Pl 1 PnJ 5 P) ( [Pl 5 PnJ i Pj 

1 ^P[l,n]-P\ [nE^^[p^,p,,])~Gn{[pi,Pn])] + 



n + I \ p 

1 / P[l,n] • P 



1 / P[l,n] • P 



n + 1 

[(n + l)F„([pi,p„]) - G„([pi,p„])] . 



^n([Pl,P»i 



?1 + 1 \ p- 

For the second line in Equation (2-3), similar calculations lead to 

-^■^«+l([Pl,Pn],P) + ^"^^ G„+i([pi,p„],p) 



1 / P[l,ri] • P 



n + 1 \ p- 



-^Pn{[Pl,Pn]) + ^"^^ Gn([pi,p„]) 



where we used 



/ ^ P[l,n] ■ P r I I n 

(P[i,«],P)^ — — , for|p|->0. 



(3-5) 



(3-6) 



(3-7) 



(3-8) 



(3-9) 
(3-10) 



Combining Equations (3-8) and (3-9), we finally obtain Equation (3-4). This concludes the proof. 
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For Pm+i, ■ • • ,Pn 0, Equation (3-2) is revised as 

^n([Pl,P„])|p^^^^...,p^^O ^ ^2 j^;„([pi,p™]), (3-11) 

where n > m. Thus, we find that Fn is approximated using Fm at a lower order m than n. The case of 
m = 1 leads to the Zel'dovich approximation in the limit that the total momentum is pi: 

F.iP.,Pn])U...,.^o ^ ^ (^) (^) • (3-12) 
Therefore, we interpret the approximation of Equation (3-2) as an extension of the Zel'dovich approximation. 



4. Reg PT 

Now, we are ready to derive Reg PT (Bernardeau et al. 2008, 2012; Taruya et al. 2012). We give an 
alternative explanation for Reg PT by the relation between SPT and the P-expansion [Equation (2-5)] and 
the approximation of the kernel functions [Equation (3-2)]. 

(r) 

Before we proceed, we show that the rth-ordcr coefficient of the P-expansion at the n-loop level P„_ioop 

(r) 

is approximately represented by PJjj.ioop at a lower m-loop level than the n-loop level: 

ri:Lp([ki,k.])|^^^^,...^^„^„ ^ 5 i^-^j p.tLp([ki,k.]), (4-1) 

where cr^ is the velocity dispersion of dark matter: cr.y = J dpPL{p)/6Tr. This is shown by substituting the 
following approximated kernel functions into Equation (2-5): 

+2„([ki,k^],pi, -pi, • • ■ ,P„,-P„)|p^^^^ ... p^^^o 

(r + 2m)! f^i^yPrn+^V f k[,,.yp„ y ffk k 1 n n n n ^ 

(r + 2n)! V P-m+i ) \ Pn J 

(4-2) 

To lead Reg PT at the one-loop level, we need the cases of m = and m = 1 m Equation (4-1): 

rilop([ki,k.])|^^...^^^^„ = -,\-^) rl:Lp([ki,k.]), 

1 / 7 2 2 \ 

ri:Lp([ki,k.])|^^...^^^^^„ = ^(-^) rW([ki,k.]). (4-3) 

From Equation (4-3), we derive an approximated expression of P|,'^ioop as follows: 

ri:Lp([ki,k,.]) ^ '^rMoop([ki,k,.])|^^^,,,^^^^„ -{n- l)r|lop([ki,k.])|^^ 

= S (-^)""r(:Lp([k.k.]) - (-^)"p<;i([k„k.]), 

(4-4) 
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where we multiply the first term in the right hand side by n because we choose a Fourier mode pi from n 



Fourier modes {i 



first term riF. 



('■) I 

n-looplP2,...,p 



1, • • ■ ,ri}. Furthermore, we need to subtract (n — l)r|i'iioop|pi,- -,p„-^o because the 
integrates the same region defined by pi , . . . , p„ — > n times in the multiple 



integral in Equation (2-5). This expression is satisfied even in the case of n = 0: F, 
Equation (4-4) into Equation (2-5), we find Reg FT with the one-loop corrections: 

F('-)(2,[ki,k,]) 



(r) 

0-loop 



Substituting 



-> 



r+2n 



n=0 



rl:Lp([ki,k,.])- 



n - 1 



exp 



F(l([ki,k,])+ (z?2rWoop[ki,k.] 



2 



r^l([ki,k,]) 



(4-5) 



As in the case of Reg FT at the one-loop level, we are further able to show the following expression: 

r|:Lp([ki,k.]) 



1 



(n-2)! 
1 



(n- 1)! 



Pal 
2 

e<7. 



r^:Lp([ki,k.]) + ( ^1 r(:Lp([ki,k.]) + ) 



rl:Lp([ki,k,.] 



1 ( k^(TlW(r) 



-Ttrccdkl, kj.]). 

This results in Reg FT with the two-loop corrections: 
FW(z,[ki,k,]) 



^)F(::i([ki,k. 



exp 



Ft;L([ki,k,]) + ( D'rtL^iikM) + :^!:^F|:l([ki,k.]) 



^^-'^^ - - - <'^) rl:Lp([k.k.]) + ) F(;i([k„k.]) 



i?^F(:Lp([ki,k,]) 



(4-6) 



(4-7) 



5. Extension of Reg PT 

In this section, we propose a natural extended version of Reg FT. In Section 4, we only applied the 
approximated kernel functions to multiple integrals involved in the coefficients of the F-expansion, such as 
the integrals of pi, . . . , p„ in Equation (2-5). Here, we further use the approximated kernel functions to 
calculate the multiple integrals in the rth-ordcr contribution to the power spectrum in the F-expansion, such 
as the integrals of ki, . . . , k^ in Equation (2-7). In other words, we apply our kernel approximation method 
to compute high-order terms of the F-expansion with Reg FT. 

The authors in Taruya et al. (2012) computed Pp'^'' up to the third order. The terms at the first, second, 
and third order in the F-expansion were calculated using Reg FT at the two-loop, one-loop (Equations (4-7) 
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and (4-5)), and tree level, respectively (see Taruya et al. 2012, Equations (23)-(26)). On the other hand, we 
evaluate any order of the F-expansion using either Equation (4-5) or Equation (4-7) to keep consisteney of 
the approximation. 

The important relation to extend Reg PT is that the rth-order coeffieient in the F-expansion is approx- 
imated by one at a lower order m than r: 



rilopdki,]^.])! 



ni+l 



'^m+1 



k[l.m] ■ k„ 



(5-1) 



where we used the following approximated kernel function, 

Fr+2n([kl,kr],Pl, -Pl, . . . ,P„, -Pn)|i, 



(m + 27i)! /k[i_„j] • km+i 



(r + 2n)! V fc^+i 



kfl.rnl ■ k„ 



-Frn+2n([kl, k„i]). 



(5-2) 



5.1. Extension of Reg PT: One-loop Level 

First, combining Equations (2-7) and (4-5), the first-order contribution to the power spectrum in the 
F-expansion with Reg PT is given by 

>n 2 



P^^\z,k)^ [F(i)(z,fc)l FL(fc) ^ 6-'='-°' 



1 + ^M-Lp(fc) 



(5-3) 



Next, we consider the second-order contribution calculated from Equation (4-5): 



2! 



(27r)3 J {27t)3 



(27r)3<5^(k-k[i,,]) [F(^)(z,ki,k2)l Pdki)PL{k2) 



d'k, rf'^2^2^^3^^(k-k[i.2])i?2PL(fcl)^?'PL(fc2) 



(27r)3 (27r)3 

(r[Z{kMy + 2F(l(ki,k2) (D'r'-;tU{k^,k,) + 'f^riZikM 



^rt^ie(kl, k2. 



(5-4) 



The first term in the square bracket leads to 



(27r)3 (27r)3 



:i27rfSD{k - k[i.2]) [F2(ki,k2)]' D^Pi^iki)D^Pi^ik2) 



^D''P22{k). (5-5) 



If we only calculate Equations (5-3) and (5-5), we obtain a result called the one- loop Reg PT solution in 
Taruya et al. (2012): 



^Rcg,l-loop(^5 fc) — 6 



-k-'D-' 



i + ^'rW,p(fc) 



D'P^{k) + e-'^-^ ""i5"P22(fc). (5-6) 



However, we further calculate the other terms in Equation (5-4). Although we need higher order terms in SPT 
than the one-loop level to strictly compute these terms, we here approximate them using P13 = 2PLri 



(1) 



■loop 



- 9 - 



and P22, which are at the one-loop leveL The second and third terms in the square bracket in Equation (5-4) 
are approximated for fci — > as 



2r^l(ki,k2)(i?^riX,p(ki,k2) + 



d^ki (Pk2 

(27r)3 (27r)3 



(2^)35c(k - ki - k2)| , _^.D^PL{kl)D^PL{k2) 



2t.(2) 



^Mlop(ki,k2) + ^^rLl(k„k2; 



<T.^2! '^^'^^ ^^^^ 



1 /ki • k2 



kl 



(27r)3 (27r)3 



2T[llik2)[D'TY'{k2) + 



(2^)3fo(k - k2)i?2PL(fcl)i?2PL(fc2) 



2t.(1) 



ri;4(fc2; 



1-loop 



(fc2) 



^ D'P^{k)-D'P^ik) 



(5-7) 



('2') l'2'l 

where we used the approximations for r^ /^^p and FJ^^q in Equation (5-1) and (5D(k— k[i^2])|fei^o = '5D(k— k2). 
Since we also need to consider the contribution from the region of fc2 ^ 0, we multiply Equation (5-7) by a 
factor of two. Equation (5-4) then becomes 



Pl^^\z,k) ^ e-'^^^^'^^i?^ (F22(fc) - fcV^PL(fc)) 



+e-— ^(fc^i?V„^)(i + i?2r(^) (fc) 



^ D'P^k). 



(5-8) 



We apply the same analysis as above to an arbitrary order coefficient of the F-expansion. The rth-order 
contribution to the power spectrum in the F-expansion is given from Equations (2-7) and (4-5) by 



P^''\z,k) e 



dPki d^kr 



(27r)3fe(k - k[i^r])D^PL{kl) ■ ■ ■ D'-PLikr) 



(27r)3 (27r)3 

F|:l([ki,k,])l' +2F|l([ki,k.]) ( Z?^Fi:/„,p([ki,k,]) + '^^^Ti:U[ki,K] 



k^D'a?. 



tree U"-! I '^r\ 



(5-9) 



The first term in the square bracket is approximated as 



d^k, d^k,. ^2^)3^^(k - k[i.,])i?2pL(fci) . • . D'P^ikr) 



rir — 1) 



r^l([ki,k.]) 



(27r)3 (27r)3 
2 



fe3,...,fer-S-0 



+ (-r(r-l)+r) F|::i([ki, k.]) 



fc2....,fcr-!-0 



(r-2)! 



-D^ (P22(fc)-fc2a2pL(fc))+e- 



2 n2^2\r-l 



(r-1)! 



-D^Pi^ik), (5-10) 
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where we multiply the first term by a factor of r(r — l)/2 because we choose the two Fourier modes ki,k2 
from the r Fourier modes, {i = 1, • • ■ ,r}. We need the second term to exclude the redundant integrated 
regions defined by k2, . . . ,kr — > 0. We used Equation (5-1) and the approximation of the delta function: 
SuCk - k[i^,r])|fc3,....fe,.^o = 5u(k - k[i^2]) and So(k - k[i^r])|fe2,...,fc.^o = fc(k - ki). The approximated 
solutions of the other terms in Equation (5-9) are also calculated by the same derivation used in the case of 



p(2). 



(r-1) 



1 + (fc) 



D^Pi^{k)-D^Pi^{k) 



(5-11) 



Therefore, we arrive at the approximated expression of Pp 



(r). 



(r-1)! 

(R^V|p^, (^-(fc) - k'alPUk)) 



D^Pi^ik) 



(5-12) 



We propose this expression as an extended version of the one-loop Reg PT solution. This certainly recovers 
the results in the cases of r = 1 and r — 2. 



Note that if we sum up all orders of the F-expansion, we derive the following power spectrum from the 
extended Reg PT at the one-loop level: 



P(z,fc) =J2Pr \^^k) = D^Pi^{k) + D^P,_,,,p{k) + (z?2fW ^p(fc) 



k^D^ar ' 



D^Pi^ik), (5-13) 



where Pi_ioop is the correction term at the one-loop level in SPT, defined as Pi_ioop = P22 + Pi3- 



5.2. Extension of Reg PT: Two-loop Level 



Similar to the derivation used in Equation (5-12), we derive the following approximated P^^^ from 
Equation (4-7): 



(r-1)! 



1 + ( D'r[]l(k) 



Ad) , . k'^D'^al o (11 , . 1 fk^D'^at 



D^PL{k) 



+,-k^D^^l {k^D^alr-^ ^, [P22(fc) - k^alP^ik)] 



+e 



(r-2)! 



(r-2)! 
(r-2)! 



D"" [P2A{k) + k'alP22{k) - Pal {P.^ik) + k^alPi^{k))] 



PiAa{k) + -^P2A{k) + ^-^P22{k) 



-k'<yt r- (fc) 



a(fc) 



+e- 



(r-3)! 



P33bik) - k^alP22{k) + &il!pL(fc) 



(5-14) 
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where P24, P336, and P^^a sltc defined as 



P2i{k) 

P44a{k) 



= 24 

= 6 
= 72 



(27r)3 (27r)3 {2nf 



(27r)-^5D(k-k[i,2])i^2(kl,k2)i^4(kl,k2,p,-p)FL(p)PL(fcl)PL(fc2), 



d^ki d^k2 d^ks 
(27r)3 (27r)3 (2^ 



■ (27r)3(5z, (k - k[i,3] ) [F3 (ki , k2 , k3 )] ' Pl (fci )Pl (fc2 )Pl (fc3 ) , 



fi3fci ^3/02 d^pi d^P2 .3„ , 



2]; 



(27r)3 (27r)3 (27r)3 (27r)3 " 

xF4(ki, k2, pi, -pi)F4(ki, k2, P2, -p2)PL{kl)PL{k2)PL{pi)PL{p2)- 



(5-15) 



Equation (5-14) recovers the two-loop Reg PT sohition which is given by (Taruya et al. 2012) 



-PRcg,2-loop(2, k) 



1 



+ pMiop(fc)+ 



+e 



2 i-ioopv / ' 2 

i^^P22(fc) + (P24(fc) + k^<jlP22{k)) 



D^Pi^ik) 



( Puaik) + -^P2A[k) + i^4^P22(fc) 



-k'D'at 



P33b{k). 



(5-16) 



In the case of ?■ = 1. Equation (5-14) coincides with the first term in Equation (5-16). For ?' = 2 and 
r = 3, we have additional terms compared to the second and third terms in Equation (5-16). This is because 
Equation (5-14) is computed using Reg PT at the two-loop level (Equation (4-7)), although the second 
and third term in Equation (5-16) are calculated by Reg PT at the one-loop level (Equation (4-5)) and 
the tree level, respectively. Furthermore, we can also derive higher order terms than the third order in the 
F-expansion from Equation (5-14). 

The power spectrum including all orders of the F-expansion is given by 

P(z,fc) = D^PL{k) + D^Pi.,,op{k) + D''P2.ioop{k) 

kW, 



8 hrW 



l-loop(^) 
P44a(fc) 



rLioopC^) 



1 fkW 
2 



P22(fc)-fcV.'(FW,pW 



^L(fc) 

D^Pl^ik) 



+D 



10 T-(l) 



2-loop 



(fc) 



fcV2(i) UkW 



PL{k), 



t2 i-ioop 2 

where P2-I00P is the correction term at the two-loop level in SPT, defined as 

P2 -loop 

= Pi5(fc) + P24(fc) + P33a(fc) + P336(fc), 



(5-17) 



(5-18) 



with Pi5 = 2F2^j^^^pPL and P33a = [F^^^iQ^pJ^pL- As in the case of the extended 1-loop Reg PT, we obtain the 
two-loop SPT solution with the additional correction terms. 
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Fig. 1. — Ratios of the predicted power spectra P to the smoothed hnear power spectrum without BAO 
Pnw presented by Eisenstein & Hu (1998) are plotted at z = 1.0. The predicted power spectra are computed 
by SPT (red sohd), Reg PT (blue sohd, Equations (5-6) and (5-16)), the extended Reg PT (red dashed, 
Equations (5-13) and (5-17)), and A^-body simulations (green symbols). 



6. Modified SPT 



We investigate the relation between the results from Reg PT (Equations (5-6) and (5-16)), the extended 
Reg PT (Equations (5-13) and (5-17)), and SPT. One might think that the resuUs of the extended RegPT 
give a better prediction for the nonlinear matter power spectra than those of SPT, because the extended 
RegPT has the correction terms as well as the SPT solution. However, this is not true. Figure 1 shows 
that the additional terms in Equations (5-13) and (5-17) hardly contribute to the predicted nonlinear power 
spectra (see the red solid and dashed lines in Figure 1). Therefore, we do not need to consider the partial 



2 

correction terms at higher order than the two-loop level in SPT, such as P^^a, ri^^ioopr2^[^^Qp, and ^'^^^ 



2-loop 

in Equation (5-17). Moreover, we clearly find that the predictions from Reg PT (blue lines) lose information 
on the nonlinear evolution of dark matter compared to the SPT predictions (red lines) due to the truncation 
of the F-expansion at second or third order. 

Now, we present the following modified SPT solution at the one- and two-loop level: 

PsPT4-ioop(^,fc) = e-'-^^^-g ^^T-fv' {D^P^ik) + (Pi3(fc) + k'alP^k))) 

r=l 

^^-.^u^.l g _ k-alP^(k)) (6-1) 



-PsPT,2-loop(^;, fc) 
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Fig. 2. — Ratios of the predicted power spectra P to the smoothed Uiicar power spectrum P^w are plotted 
at the redshifts of z = 3.0 1.0, and 0.35. The predicted power spectra are computed by the one-loop 
SPT (purple solid), the two-loop SPT (red solid), the two-loop Reg PT in Equation (5-16) (blue solid), 
the closure theory (orange solid), the A^-body simulations (green symbols), and the modified two- loop SPT 
(Equation (6-2)) where the finite truncated orders of the F-cxpansion go from 1 to 8 (black dashed; 7'^ = 1-5, 
7, and 8 from left to right in order; black solid: =6). Furthermore, the fractional differences, defined 
as Diff[%] = (PNbody — P) * 100/PNbody where PNbody is the predicted power spectrum by the TV-body 
simulations, are plotted at each redshift. There, the plotted models are the one- loop SPT (purple symbols), 
the two- loop SPT (red symbols), the two- loop Reg PT (blue symbols), the closure theory (orange symbols), 
and the modified SPT truncated at the sixth order of the F-expansion (black symbols). 



-k^D 



/'I,2 n2„2^r-l 



D'P^ik) + (Pi3(fc) + k^alPLik)) 



E 

r=2 



(r-2)! 



+ 0" P33a{k) + Pl^ik) + fcV^Pi3(fc) + ^—^P^{k) 



(P22(fc) - ealPi^{k)) 



D"" {P2i{k) + k'alP22 - k'al {P,3{k) + k^alPi^ik))) 



,,-.wg- (w^gr-^ ^ar„ .2.2„ , (fcv„^)^ 



r=3 



(r-3)! 



PMk) - fc^<P22(fc) -t- ^^r^PLik) 



(6-2) 



These expressions are obtained by removing the additional terms at higher orders than the one- and two-loop 
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level in SPT from Equations (5-12) and (5-14), respectively, and cutting off the F-expansion at a finite order 
Tfc. For rfc oo, the above solutions reduce to the usual SPT solutions. Note that these solutions avoid the 
divergent behavior of the power spectrum at small scales by considering the finite order one needs. 

In Figure 2, we show a comparison between the predicted power spectra from the modified SPT solution 
at the two-loop level (Equation (6-2)) with the finite truncated orders of the F-expansion going from 1 to 8 
(black dashed), the one- loop SPT (purple solid), the two- loop SPT (red solid), the two-loop Reg FT (blue 
solid), the closure theory (Taruya et al. 2009) (orange solid), and the A^-body simulations (green symbols). 
First, the predictions of the two-loop SPT agree better with the A^-body simulations than those of the 
two-loop Reg FT and the closure theory at the redshifts {z = 3.0, 1.0, and 0.35). Next, the modified SPT 
solutions truncated at the third order of the F-expansion give almost the same results as the two-loop Reg 
FT. As previously mentioned, this is because the additional correction terms in the two-loop Reg FT have 
no effective information on the nonlinear evolution of dark matter. Third, the divergent behavior in the 
two-loop SPT at small scales is indeed modified by the truncation of the F-expansion (see the black dashed 
lines). We emphasized the sixth-order solutions of the F-expansion in the modified SPT using the black 
solid lines. The sixth-order solutions accidentally agree with the A^-body results better than the two-loop 
SPT solutions due to the damping factor from the truncated F-expansion (see the red and black symbols in 
the fractional differences in Figure 2). We suggest that this solution should be used instead of the two- loop 
SPT. The solution guarantees accuracy comparable to the two-loop SPT and has no divergent behavior at 
small scales. 



7. Conclusion 

Based on the idea that any resummation theory should be described and solved through the kernel 
function in SPT, we established an approximation of the kernel function (Equation (3-2)). 

This approximation explains the existing Reg FT model, and gives its natural extension (Equation (5- 
14)). However, the extended Reg FT is equivalent to SPT with negligibly small correction terms, and gives 
almost the same result as for the SPT (see Figure 1). Therefore, we do not need to consider these additional 
terms. Moreover, since the Reg FT solutions (Equations (5-6) and (5-16)) truncate the F-expansion at the 
second and third order, they do not contain more effective information on the nonlinear dark matter evolution 
than the SPT solutions. In other words, while the original Reg FT solution becomes zero beyond BAO scales 
due to exponential damping, the extended Reg FT solution allows us to predict the power spectrum at the 
higli-fc region, because the exponential damping factor is completely canceled out. In fact, the nonlinear 
power spectra computed in the two-loop SPT agree better with the A^-body results than that of the two-loop 
Reg FT at redshifts of z = 3.0, 1.0, and 0.35, especially at high-fc region. Hence, the two-loop SPT is one of 
the best models to predict the nonlinear matter power spectrum even at present (see Figure 2). 

In the spirit of the F-expansion, each loop corrections are concentrated at a narrow range of scales, 
and the reliable range of the approximation becomes much more obvious. Therefore, the Reg FT two-loop 
solution, which truncates the F-expansion at the third order, guarantees the validity of the predicted power 
spectrum at small scale regions. In this paper, we showed that the Reg FT one- and two-loop solutions 
are a part of the SPT one- and two- loop solutions plus negligible terms (Equation (5-17)) and extract the 
information at small scale regions in the SPT solutions. Thereby, the predicted power spectra in the Reg FT 
model have the same behavior as those in the SPT model at the BAO scales (see Figure 2). Furthermore, 
we interpret the truncation of the F-expansion as the loss of the information at small scale regions. In this 
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sense the Reg PT solution docs not have more effective information on the nonhnear power spectrum than 
the SPT solution. 

We presented modified versions of the SPT solutions at the one- and two-loop level (Equations (6-1) 
and (6-2)), which are the SPT solutions described in the context of the F-expansion. These solutions avoid 
the divergent behavior of the SPT solutions at small scales by the truncation of the F-expansion at the 
finite order one needs. We suggest that one should use the modified SPT solution at the sixth order of the 
F-expansion instead of the two-loop SPT solution, because this solution guarantees accuracy comparable to 
the two-loop SPT and has no divergent behavior. 

Finally, we shall mention some applications of our result. So far the correlation function has not been 
computed from SPT because of the divergent behavior. However, the modified SPT solution would allow 
us to do it. Similarly, we occasionally encounter the difficulty that we have to compute the integral of the 
linear or nonlinear power spectrum but the integrand does not converge, such as the calculation of back 
reaction from the matter perturbations. Our result would be useful for this case. However, we leave these 
applications as our future work. 
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